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A global monopole (or other topological defect) formed during a recent phase transition with core 
size comparable to the present Hubble scale, could induce the observed accelerating expansion of the 
universe. In such a model, topological considerations trap the scalar field close to a local maximum of 
its potential in a cosmologically large region of space. We perform detailed numerical simulations of 
such an mhomogeneous dark energy system {topological quintessence) minimally coupled to gravity, 
in a flat background of initially homogeneous matter. We find that when the energy density of 
the field in the monopole core starts dominating the background density, the spacetime in the core 
starts to accelerate its expansion in accordance to a ACDM model with an effective inhomogeneous 
spherical dark energy density parameter Q\(r). The matter density profile is found to respond to the 
global monopole profile via an anti-correlation (matter underdensity in the monopole core). Away 
from the monopole core, the spacetime is effectively Einstein-deSitter (Q.A{r out ) — > 0) while at the 
center Oa^ — 0) is maximum. We fit the numerically obtained expansion rate at the monopole core 
to the Union2 data and show that the quality of fit is almost identical to that of ACDM. Finally, 
we discuss potential observational signatures of this class of inhomogeneous dark energy models. 

PACS numbers: 98.80.Es,98.65.Dx,98.62.Sb 



1. INTRODUCTION 

The simplest cosmological model consistent with the 
observed accelerating expansion of the universe is the 
ACDM model [1] which assumes that about 70% of the 
energy density of the universe consists of a homogeneous 
and constant in time energy density the cosmological con- 
stant. Despite of its simplicity, this model is faced by the- 
oretical (fine tuning and coincidence problem) and obser- 
vational challenges (eg the CMB anomalies [3] and large 
scale velocity flows [4, 5]). 

Most models generalizing ACDM are motivated by the 
theoretical challenges of the model and they are based 
on breaking the time translation invariance of the cos- 
mological constant while retaining large scale homogene- 
ity. They involve either dark energy [1, 2] (a time de- 
pendent form of energy with negative pressure) or ex- 
tensions of general relativity with repulsive gravitational 
properties on large scales. However, observational chal- 
lenges of ACDM (CMB anomalies and large scale veloc- 
ity flows) appear to hint towards violation of large scale 
homogeneity [6]. Theoretical models generalizing ACDM 
and involving violation of large scale homogeneity include 
Lemaitre-Tolman-Bondi (LTB) models [7] based on large 
voids of matter and no dark energy [8] , anisotropic dark 
energy [9, 10], fundamental modifications of large scale 
topology and geometry [11], a Hubble scale magnetic field 
[12], dark energy with Hubble scale inhomogeneities [13] 



'Electronic address: jcbuenoSfis.ucm.es 
^Electronic address: leandrosSuoi.gr 



etc. A common problem of these models is that most of 
them lack a well defined causal mechanism that gener- 
ates the proposed Hubble scale inhomogeneities starting 
from early time homogeneity. For example the size of the 
matter void required by LTB models is about 2Gpc and 
is unlikely to produce in the context of any known viable 
cosmological model. 

In this study we extend recent work investigating 
observational constraints on the existence of a Hubble 
scale spherically symmetric dark energy inhomogeneity 
[13]. Such inhomogeneity was found to be consistent 
with current observations provided that its radius is 
larger than about 2Gpc [13]. We point out that such 
inhomogeneity can be naturally produced by a global 
monopole with a Hubble scale core formed during a re- 
cent phase transition. Such a mechanism, called topologi- 
cal quintessence has similar features as 'topological infla- 
tion' [14-16] which may have occurred in the early uni- 
verse and involves accelerating cosmological expansion 
triggered in the core of global topological defects when 
proper conditions are satisfied. The investigation of these 
conditions and the required values of model parameters 
in the context of the recent universe where matter is also 
present is the goal of the present study. 

The structure of this paper is the following: In the 
next section we set up the basic equations of gravitating 
global monopoles in an expanding universe with mat- 
ter and use analytical approximations to derive the pa- 
rameter values that would lead to accelerating expan- 
sion in the monopole core at recent times and on Hubble 
scales. In section 3 we present a numerical solution of 
the full system of field equations coupled to gravity and 
show that indeed accelerating expansion takes place in 
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the monopole core when its energy density starts domi- 
nating over matter. We also fit the predicted expansion 
rate in the core to the Union2 data. Finally in Sec. 4 we 
summarize our main results and discuss future prospects 
of this project. 



GRAVITATING GLOBAL MONOPOLES AND 
COSMOLOGICAL EXPANSION 



the model r\ (the vacuum expectation value) and A (the 
coupling constant). The global monopole field configura- 
tion is described by the hedgehog ansatz 
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<f> a = $(r, t)f a = $(r, t) (sin cos 93, sinflsintp, cos( 
with boundary conditions 

$(0,t) = , $(00,*)=^, (2.4) 



Global monopoles [17] are formed in the context of 
scalar field theories when the vacuum manifold M. con- 
tains surfaces that cannot be continuously shrunk to 
a point namely when tt 2 {M) 7^ The correspond- 
ing global field configurations are topologically distinct 
from the vacuum and have most of their energy localized 
in a spherically symmetric region of space, the global 
monopole core. The vacuum energy localized in the 
global monopole core has negative pressure and there- 
fore can induce accelerating cosmological expansion. In 
order for this expansion to have interesting observable 
cosmological consequences two basic conditions must be 
satisfied: 

• The scale of the monopole core should be cosmo- 
logical and comparable to the Hubble scale. This 
condition is necessary so that the effects of the ac- 
celerating expansion are applicable on cosmological 
scales thus affecting cosmological observations (eg. 
Hubble diagrams). 

• The vacuum energy in the core dominates over 
other forms of energy with non-negative pressure 
(matter or radiation). This condition is necessary 
so that cosmological dynamics is determined by the 
negative pressure of the monopole rather than by 
the other forms of energy present. 

These two conditions are easily translated to constraints 
on the parameters of the field theory that predicts the 
existence of the global monopole. Consider for example 
the action 
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where C m is the Lagrangian density of matter fields, 
$ a (a = 1, 2, 3) is an 0(3) symmetric scalar field and 

1 



y($) = -A($ 2 -t] 2 ) 2 , $ = V$ a $ a . (2.2) 

The vacuum energy density in the monopole core and the 
size of the core are determined by the two parameters of 



where we have allowed for a time dependence having in 
mind a cosmological setup of an expanding background. 
The general spherically symmetric spacetime around a 
global monopole may be described by a metric of the 
form 

ds 2 = -dt 2 + A 2 (r, t)dr 2 + B 2 (r, t)r 2 (d9 2 + sin 2 Odcp 2 ). 

(2.5) 

In order to derive the cosmological dynamical equations 
we need the total energy momentum tensor T av . We 
assume a cosmological setup at recent cosmological times 
and therefore T^ v is dominated by the global monopole 
vacuum energy and matter, ie. 



rp rp(mon) 1 rp(mat) 

1 V-V ~ 1 UV T 1 III/ 
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The energy momentum tensor of the monopole is given 

by 
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For example for the monopole energy density we have 

£2 A , 2 _ 2V 

2 + 2A 2 + B 2 r 2 + 4 ^ V > 

(2-8) 

The energy-momentum tensor of matter is that of a per- 
fect fluid with zero pressure (p = 0). 
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where p mat 
the velocity 4-vector of the fluid 

= ( 



,0,0). (2.10) 



with v = v(r, t) the radial matter fluid velocity. Extrem- 
ising the action (2.1) with the metric (2.5), the ansatz 
(2.3) and the presence of matter we obtain the dynami- 
cal equations 
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where prime (dot) denotes differentiation with respect to 
radius r (time t) and 
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K = Xj. (2.15) 



The conservation of the energy momentum tensor for the 
matter fluid (2.9) in the metric (2.5) leads to the following 
equations for v and p mat 
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Notice that in the homogeneous FRW limit (B(r,t) — 
A(r,t) = a(t), $>(r,t) = rj) we obtain the familiar Fried- 
man equations with v — and the usual conservation of 
matter equation. 

It is straightforward to show [17] that for p mat = 0, the 
system (2.11)-(2.14) has a static solution <&(r, t) = f(r) 
that is stable for a range of the parameters A and rj of the 
model. It is straightforward to show using (2.14) that in 
a locally flat space (A ~ B ~ 1), the core scale (where 

$/?7 ~ 1) of this static global monopole is 



5~ A 
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while the vacuum energy density in this core region is 
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In addition, in the flat space approximation the asymp- 
totic behavior of the static global monopole solution may 



be found using Eq. (2.14) to be of the form 
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r»J. (2.21) 



The above described flat space approximation is valid as 
long as rj <C mpi . For rj comparable to mpi the effects of 
gravity must be taken into account and the coupled sys- 
tem (2.11)-(2.14) needs to be solved. An approximate so- 
lution may be obtained valid far away from the monopole 

2 

core [18] (r 3> 8) where <&(r) ~ r\ and p mon ~ 2 . In 
the coordinates of the metric (2.5), this approximate so- 
lution is of the form [18]: 

(2.22) 
(2.23) 

where the contribution of matter has been ignored. In the 
Appendix we show a generalization of eqs (2.22), (2.23) 
when matter is also present. 

Deep inside the monopole core (r < i), spacetime is 
approximately homogeneous and we may set A(r, t) ~ 
B(r,t) = a,i n (t). In this limit, Eq. (2.11) becomes 
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where ao,i n = din(to) is the scale factor at the present 
time to, P§in is the present density of matter in the 
monopole core and 



H in (t) = ( ^ 

. Q>in 



(2.25) 
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Similarly, far away from the monopole core (r 3> 5), 
spacetime is also approximately homogeneous and we 
may set A(r, t) ~ B(r,t) = a out (t). In this limit, 
Eq. (2. 1 1) becomes 
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corresponding to a flat matter dominated universe. 

In the context of topological quintessence a cosmolog- 
ically large region of space associated with the core of 
a recently formed global monopole starts an accelerat- 
ing expansion in accordance with Eq. (2.24). For proper 
values of the parameters of the model A and r\ (so that 

~~4~ ~ Po'tn) the local expansion rate predicted by (2.24) 
is identical to that of ACDM. 

Around a time t when Ar/ 4 3> p™ at the vacuum energy 
dominates over matter, the scale factor may be approxi- 
mated as 



a(t) ~ a(t n )e Ha{t '^ 



where 



H = H(to) 
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The slow evolution of the scalar field $ may also be ap- 
proximately obtained using (2.14). For r -c 6 and in the 
slow roll approximation this becomes 



3JT(t )$ A?7 2 $, 
where t is the present time. This leads to [14] 

$(r,t) ~ $(r,t )e^ ( * _to) . 



In the limit 
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(2.30) 



(2.31) 



the field is practically static on cosmological timescales 
(At ~ Hq 1 ) and the static monopole approximation of 
Eqs. (2.18)-(2.21) holds. Notice however, that despite 
of the approximately static profile of the scalar field, 
the expansion rate inside the monopole core (Eq. (2.24)) 
may differ significantly from the expansion rate far away 
(Eq. (2.26)). 
For 
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the field profile evolution is significant on cosmological 
timescales. The corresponding condition on the parame- 
ters of the model is easily obtained from Eqs. (2.31) and 
(2.28) as 



mpi ~ 0.4 mpi . 



(2.33) 



This is the condition for breakdown of the static 
monopole solution. A similar result (derived in a dif- 
ferent manner) was also presented in Ref. [14] in the 
context of topological inflation (see also [15]) and con- 
firmed numerically in Ref. [19] and with the presence of 
radiation in Ref. [18]. The numerically obtained condi- 
tion for accelerating expansion of the monopole core (in 
physical coordinates) was found to be n > 0.33 mpi in 
good agreement with Eq. (2.33). Notice that in this case 
there is no condition imposed on the coupling constant 
A. 

In the context of late time accelerating expansion in the 
presence of matter, the two conditions that need to be 
satisfied for consistency with observations are translated 
to conditions on the model parameters as follows: 

1. The scale of the monopole core should be compa- 
rable to the Hubble scale (5 — aH^ 1 where a is a 
parameter whose range is to be determined obser- 
vationally). This is expressed by the condition 
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where Oom = ^ — is the present fractional den- 
sity of matter over the critical density p 0cr ~ 
8 x 10~ 47 h 2 GeV 4 inside the core, with h the Hub- 
ble parameter in units of 100km/ (sec ■ Mpc). 

2. The vacuum energy density in the monopole core 
must be large enough to start dominating over the 
matter energy density today. This condition im- 
plies that 

A ?y 4 

— ^— = &,TQ POcr = (1 — ^Om) POcr 

~ 4 (1 - f2 0m ) 10~ 123 m 4 Pl h 2 , (2.35) 

where £Itq = 1 — ^om — 0.73 is the normalized 
vacuum energy density in the monopole core. 



Using Eqs. (2.34) and (2.35) we find 
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(2.36) 
(2.37) 



where fj = is the dimensionless vacuum expectation 
value of the global monopole field potential. The fine 
tuning of the coupling constant A is worth noting but it 
is not worse than the fine tuning corresponding to the 
standard ACDM model. 

We have therefore used analytical arguments to derive 
the values of the global monopole model parameters that 
are required to induce the observed accelerated cosmolog- 
ical expansion in the context of topological quintessence. 
In what follows we solve the system (2.11) - (2.17) numer- 
ically in order to demonstrate explicitly that the topolog- 
ical quintessence mechanism is observationally viable. 
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3. NUMERICAL SOLUTION OF MATTER - 
GLOBAL MONOPOLE SYSTEM 

A. Method 

In order to solve the system (2.11) - (2.14) numerically 
we use the following rescaled dimensionless variables: 



r = VXr/r , t = V\rjt 



-mat 



At? 4 



-mon — P 



At, 4 



and in what follows we omit the bar for simplicity. It 
is straightforward to express the system of Eqs. (2.11)- 
(2.17) in terms of the above dimensionless quantities. 
The form of the equations remains almost unchanged 
with the following minor modifications: 

• We replace all dimensional parameters with the cor- 
responding barred dimensionless while A and 77 get 
eliminated by the rescaling and are replaced by 1 
in all equations. 

• The factor -^5- on the righthand side of Eq. (2.11) 

gets replaced by 87T7/ 2 , and similarly for Eqs. (2.12) 
and (2.13). 

We solve the system of Eqs. (2.12)-(2.17) (in its rescaled 
form) and use the constraint in Eq. (2.11) as a test of the 
accuracy of the solution. We find that in all runs this 
constraint is satisfied at a level better than 1%. We start 
the evolution at a rescaled time to well in the matter era 
assuming an initially homogeneous background of matter 
p mat (r,t ) = 4.5, which implies that 



'Mo) 



= 18 



(3.1) 



p ^.<. = p "»<"»(o,t) = I. The initial field profile 
corresponds to the static monopolc solution <&(r, to) = 
f(r). We thus use the following initial conditions chosen 
so that Eq. (2.11) is satisfied (all quantities are rescaled 
dimensionless): 



P mat (r,t ) 
v(r, t ) 



to = 
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Y 6p(t ,r)Trr] 2 
B(r,io) = l, 
f(r), 
0, 



(3.2) 
(3.3) 

(3.4) 

(3.5) 
(3.6) 
(3.7) 



- Vp mQt Mo ) + P mon Mo) , (3. 
37? 



K{r,t ) = p mQt Mo) + P mon {r,to) . (3.9) 



These initial conditions correspond to the physical as- 
sumption of the formation of a global monopole during a 
recent phase transition that takes place in a fiat homoge- 
neous matter dominated universe. Initial curvature could 
be introduced by modifying the initial condition for the 
scale factors A(r,to) = B(r,to) = 1. The corresponding 
boundary conditions used during the evolution are the 
following: 



$(0,t) 
*'(0,t) 
K(0,t) 



0, 

f(r = 0) , 
3 iff (0,t). 



(3.10) 
(3.11) 
(3.12) 



The present time t p where the evolution of the above 
system stops is determined by demanding that 
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0.37. 
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where we have set Slom — 0.27. It is now straightforward 
to evolve the above cosmological system from the initial 
time to to the present time t p . In the next subsection we 
present the results of this evolution and the comparison 
with observations. The key assumption in the interpre- 
tation of these results is that the core size of the global 
monopole constitutes a large part of the present Hubble 
scale as discussed in the previous section. 



B. Results 

The only parameter that needs to be fixed in the evo- 
lution of above system is the scale of symmetry break- 
ing 77. We have considered values of 77 in the range 
77 G [7 x 10~ 4 , 0.6] with no significant change in the form 
of the rescaled expansion rate profile at recent cosmo- 
logical times. The results presented here correspond to 
77 = 0.1 unless otherwise noted. The value of the coupling 
constant A is fixed implicitly for each 77 by demanding 
that the energy density at the monopole core is compa- 
rable to the present matter density 



p mat (0,t p ) = 0.37 p core = 



0.37 



(3.14) 



(see Eq. (3.13)). 

In Fig. 1 we show the evolution of the scalar field profile 
for 77 = 0.1 and 77 = 0.6 from the time to up to the present 
time t p . As expected from Eq. (2.33), for 77 = 0.1 the 
field rolls down to its vacuum on a timescale less than 
Hq 1 . This roll-down takes place at all comoving radii r 
and manifests itself as a collapse of the scalar field profile 
in comoving coordinates. For 77 = 0.6 the effect is much 
slower because for constant p core (as is the requirement 
is our case) the exponent of Eq. (2.30) is written as 



At? 2 
3ffo 



At 



4 p c 



At 



3 H if 



(3.15) 



since in our case At = t p — to, Ho and p core are fixed 
as we vary 77. Thus as we increase 77 by a factor of 6, 
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the exponent of Eq. (2.30) decreases by a factor of 36 
and the evolution of $ gets suppressed. Indeed we have 
checked that in our simulations the rescaled exponent of 
eq. (2.30) tp H to varies with rj as rf 2 as expected from 



the above arguments. 



by many orders of magnitude below rj = 0.1 since this 
would lead to a monopole core scale that is much smaller 
than the Hubble scale and therefore would have minor 
cosmological implications. Thus it is hard to avoid the 
fine tuning implied by Eqs. (2.36) and (2.37). 
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FIG. 1: Evolution of the scalar field profile for r\ = 0.1 (top 
panel) and rj — 0.6 (bottom panel) from the time fo up to the 
present time t p . The profiles shown correspond to to (dotted 
line), tp/2 (dot-dashed line), 2t p /3 (dashed line), t p (solid 
line) and the profile collapses as time increases. Notice that 
the evolution gets much slower as we increase rj because the 
exponent of Eq. (2.30) decreases as rj~ 2 for fixed monopole 
core density. 

In Fig. 2 we show the evolution of the scale factors 
A(r, t) (upper panel) and B(r 7 t) (lower panel) profile for 
rj = 0.1 and rj = 0.6 from the time to up to the present 
time t p . Even though the evolution of the scalar field 
changes significantly as we change rj from rj = 0.1 to 
rj = 0.6 (see Fig. 1) the corresponding evolution of the 
scale factors around the present time t p is mildly mod- 
ified. This is due to the fact that as we decrease rj, A 
has to increase more rapidly so that the monopole core 
contribution to the energy density of the universe at the 
present time t p remains unchanged and comparable to 
the contribution of the matter density. This rapid in- 
crease of A leads to a decrease of the monopole core scale 
compared to the Hubble scale. Thus, even though the 
expansion rate in the monopole core at the present time 
is not affected by the value of rj, we cannot decrease rj 
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FIG. 2: Evolution of the scale factors A(r, t) and B(r, t) profile 
for rj = 0.1 (black lines) and rj = 0.6 (blue lines) from the time 
to up to the present time t p . A mesh of dots is superimposed 
to the curve corresponding to rj = 0.1, t = t p to distinguish 
it from the curve obtained for rj = 0.6, t = t p . The profiles 
shown correspond to to (dotted line), t p /3 (dot-dashed line), 
2t p /3 (dashed line), t p (solid line). Higher curves correspond 
to more recent times. Notice that even though the evolution 
of the scalar field changes significantly as we change rj from 
rj — 0.1 to rj = 0.6 (see Fig. 1) the corresponding evolution 
of the scale factors is only slightly changed. The (rescaled) 
region of accelerated expansion is slightly smaller in the case 
of smaller rj as expected due to the more rapidly shrinking 
core seen in Fig. 1. 

It is now straightforward to obtain the Hubble expan- 
sion rate corresponding to each scale factor as a func- 
tion of redshift and of the comoving distance from the 
monopole center. For example in the case of the scale 
factor A(r, t) we have 



H A (r,z) 



A(r,t(r,z)) 
A{r,t(r, z)) 



(3.16) 



where the function i(r, z) is obtained by solving numeri- 
cally the equation A(r, t) = — j-^ with respect to t. 



7 



The derived forms of the Hubble expansion rates 
Hji(r,z), Hb(t,z) are well fit by the ACDM form 
of H A (r, zf = Hi (fi A (r) + (1 - A (r))(l + zf) where 
17a (r) interpolates smoothly between Oa(0) = 0.73 and 
^aO" 3> 1) — 0. This is demonstrated in Fig. 3 where 
we show HA{r 1 z)/HA{r 1 0) and Hs{r 1 z)j Hb{t,0) for 
r = 0, r = 0.5 and r = 5 along with the correspond- 
ing H A (r,z)/H for 17a = 0.73 (thick, dashed line), for 
Oa = (dot-dashed line) and for the best fit J7a corre- 
sponding to r = 0.5. 
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z 

FIG. 3: Ratios H A (r, z)/H A (r, 0) and H B (r, z)/H B {r, 0) for 
r = (lower solid line), r = 0.5 and r = 5 (upper solid lines) 
along with the corresponding Ha(v, z)/Hq for J1a = 0.73 
(thick, dashed line) and for Qa = (dot-dashed line). Also 
included for comparison H\(r, z)/Ho for the best fit values 
of OA(r) shown in Fig. 4: ^(0.5) ~ 0.27 (top panel) and 
S7a(0.5) ~ 0.54 (bottom panel). We have set n = 0.1. 

The profiles of f2 A (r) that provide the best fit to 
Ha{t, z)/Ha{t, 0) and to Hb{t, z)/Hb{t 1 0) are shown in 
Fig. 4. A similar profile for H B {r,z)/H B (r 7 0) has been 
used in Ref. [13] to derive detailed observational con- 
straints for inhomogeneous dark energy models. It was 
argued using qualitative arguments that the expansion 
rate predicted by topological quintessence could be well 
described by such a profile. Our numerical analysis con- 
firms this expectation and also derives the precise form 
of this profile shown in Fig. 4. 

The isocurvature nature of the global monopole dark 
energy inhomogeneity is demonstrated in Fig. 5 where 
we show the evolution of the dark energy and matter 
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FIG. 4: Values of Qa(t) that provide the best fit to 
H A (r,z)/H A (r,0) (solid line) and to H B (r, z)/ H B (r,0) 
(dashed line). 



density profiles. The dark energy density profile is de- 
fined in Eq. (2.8). Notice the anticorrelation between 
dark matter and dark energy inhomogencitics which is 
expressed by the development of matter underdensity in 
the monopole core where dark energy is overdense. This 
nonlinear result confirms the perturbative analysis of Ref. 
[20] where this anticorrelation was also pointed out. The 
evolution of the velocity of matter was also evaluated dur- 
ing the simulation and was found to remain very close to 
throughout the evolution. 

Using the derived forms of the Hubble expansion rates 
it is now straightforward to obtain the luminosity dis- 
tance using Hb{t,z). This luminosity distance may be 
fit to the Union2 data [21] and the value of x 2 ( r )/d.o.f. 
(x 2 per degree of freedom) may be evaluated as a func- 
tion of r. The predicted distance moduli curves (resid- 
ual with respect to ACDM best fit) for various values 
of r are shown in Fig. 6 (top panel) superposed on the 
Union2 data. Notice how the agreement deteriorates as 
r increases away from the core r ~ 1 to the Einstein-de 
Sitter (EdS) matter dominated region. This is quantita- 
tively expressed in Fig. 6 (bottom panel) where we show 
the x 2 { r )/d.o.f. which changes from a value less than 1 
close to the center to a value close to 2 away from the 
core. Thus a global monopole with a Hubble scale core 
provides a good fit to the Union2 data. A corresponding 
analysis of spherical inhomogeneous dark energy, based 
on a toy model but taking into account the detailed in- 
homogeneous structure of the dark energy profile may be 
found in Ref. [13]. 



4. CONCLUSION 

We have introduced a simple new mechanism called 
topological quintessence that constitutes a generic gener- 
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FIG. 5: Evolution of the dark energy (top panel) and properly 
normalized matter (bottom panel) density profiles. The pro- 
files shown correspond to to (dotted line), t p /3 (dot-dashed 
line), 2t p /3 (dashed line), t p (solid line). Notice the matter 
underdensity that develops in the monopole core while the 
monopole density appears to slowly collapse to the center in 
comoving coordinates. 



alization of ACDM. Instead of breaking time translation 
invariance of the cosmological constant which occurs in 
most studied generalizations of ACDM (eg. the usual 
quintessence) , topological quintessence involves breaking 
space translation invariance of the cosmological constant. 
This breaking is achieved via a recent low energy phase 
transition which gives rise to topological defects with 
Hubble scale core and energy density in the core com- 
parable to the present matter density. 

We have studied in some detail the case of spherical 
global defects (global monopoles) while our mechanism 
can be easily generalized to different symmetries and dif- 
ferent field varieties. We have shown that for proper- 
values of field theory parameters shown in Eqs. (2.36) 
and (2.37) topological quintessence is consistent with the 
Union2 Type la supernovae data and predicts recent ac- 
celerating expansion on Hubble scales very similar to the 
one predicted by ACDM. Due to the cosmologically large 
scale of the dark energy inhomogeneity there is also no 
breaking of the Copernican principle. 

The recent nature of the phase transition involved in 
topological quintessence allows for a matter era that is 
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FIG. 6: Top Panel: Predicted distance moduli curves (resid- 
ual with respect to ACDM best fit) based on the expansion 
rate for various values of r are shown in the left panel su- 
perposed on the Union2 data. The distance moduli curves 
correspond to r = (best fit to Union2; red line), r = 0.25, 
r = 0.5, r — 5 (worst fit corresponding to EdS matter domi- 
nated universe away from the monopole core). Bottom Panel: 
The quality of fit of the expansion rate Hb (r, z) to the Union2 
data expressed through \ 2 { r )/d.o.f.. The dashed line corre- 
sponds to ACDM in both panels. 



for the most part indistinguishable from that of the stan- 
dard ACDM model. Therefore only the large scale low 
multipole part of the CMB spectrum is expected to be 
modified. A general study of such anticipated observa- 
tional consequences in this class of models may be found 
in Ref. [13]. In order to implement the general observa- 
tional constraints obtained in Ref. [13] to the particular 
case of topological quintessence, we need to establish the 
correspondence between the general parameters of Ref. 
[13] and the parameters of the topological quintessence 
models. This correspondence may be summarized as fol- 
lows: 

• The size of the dark energy inhomogeneity ro of 
Ref. [13] should ne identified with the scale of the 
global monopole core S ~ A -1 ' 2 ^ -1 . 

• The normalized dark energy density in the center 
of the inhomogeneity (ftxin of Ref. [13]) should be 
identified with the normalized potential energy of 



9 



core \ 4 

the global monopole center — p — ~ j 2 - where p c is 
the present critical density for flatness. 

• The distance r a i, s of the observer from the center 
has an identical meaning in both approaches. 

The constraints on the parameters ro, flxin and r b s 
obtained in Rcf. [13] may be summarized as follows: 
The Union2 Snla data constrain r , Qxin = 1 — fiom 
to 15Gpc > ro > 1.76Gpc (the upper bound corresponds 
to homogeneous ACDM ) 0.8 > Qxin > 0.54 at the 3a 
level (Fig. 2ofRef. [13]). The 3(7 range of r corresponds 
to a range of the parameter a (defined above eq. (2.34)) 
as 5 > a > 0.6. Using the 3cr range of the parameters 
Qxin and r b s (or equivalcntly f^om and a defined above 
cq. (2.34)) and equations (2.36) and (2.37) we obtain 
the observationally determined 3cr range of the model 
parameter A and f). We find at the 3cx level 

4 x 10~ 117 V > A > 1.1 x 1Q- I22 h 2 (4.1) 

and 

0.63 > fj > 0.02 (4.2) 

The value of r i, s is very weakly constrained by the 
Snla data and it can be close to the tq (Fig. 6 of Rcf. 
[13]) especially for large inhomogeneity scales ro (AGpc 
or larger). The value of r bs is more severely constrained 
by the smallness of the observed CMB dipole. For r up 
to 7Gpc, r bs is constrained to be less than llOMpc (Fig. 
7 of Ref. [13]) while the CMB dipole constraint on r a i, s 
is alleviated if the inhomogeneity reaches the distance to 
the last scattering surface (ro ~ 13Gpc). An extended 
study of these constraints as well as derivation of spe- 
cific observational signatures of topological quintessence 
is an important next step but is beyond the scope of the 
present study. 

The following directions could lead to useful extensions 
of the present project: 

• A more detailed investigation of the above men- 
tioned observational consequences and their rela- 
tion to CMB anomalies [3] and other puzzles [22] 
of the standard ACDM model [1]. 

I 



• The consideration of multi-defect configurations 
could reveal interesting new effects and observa- 
tional signatures. Such signatures are expected to 
include Hubble scale velocity flows experienced by 
off center observers, aligned low multipole CMB 
moments in multi-defect configurations, directional 
variations of fundamental constants due to possible 
coupling to the defect scalar field etc. 

• The application of our analysis to the case of non- 
minimally coupled scalar field (extended topolog- 
ical quintessence). In this case we anticipate a 
more drastic reaction of the matter profile at the 
monopole core. 

• The consideration of alternative topological defects 
(a global string, a wall or a gauged defect). The 
identification of distinct observational signatures 
for each defect geometry constitutes an useful probe 
of this class of models. 

• Finally our numerical tools may also be used with 
proper initial conditions to investigate the effect 
of dark energy (minimally and non-minimally cou- 
pled) on the profiles of bound matter structures 
like clusters of galaxies. This analysis is currently 
in progress. 



Numerical Analysis Files 

The Mathematica files used to produce the 
figures of this study may be downloaded from 
http://leandros.physics.uoi.gr/topquint . 



Appendix 

In the presence of matter and far away from the 
monopole core (r 3> 5), using Eq. (2.11) we find that the 
approximate form of the scale factors A(r,t) and B(r,t) 
satisfy the equations: 



* \ y Omat y £>mat ]j Ornat J ^mat 

A(r,t) ~ B(r,t) (4.4) 

I 

where (3 = ^ eq , B eg (r, t eq ) is the scale factor at the time of the monopole (p mat (t eq ) = p mon (t eq ), B(r, t eq ) = B eq ) 
when the density of matter is equal to the energy density and B mat is the scale factor at an early time t mat when 
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matter dominates (t ma t <C t eq ). It is easy to check that 
the above expression reduces to the anticipated forms in 
the matter dominated (B <C B eq ) and in the monopole 
dominated (B > B eq ) eras. These expressions are of the 
form 

A(r,t) ~ B(r,t)~ ' , B <c B eq (4.5) 

4(r,t) ~ B (r,t)~J^-- , 5»B e9 (4.6) 
V 3 mi,, r 



where the rescaling of the scale factor in the monopole 
era ^= compared to Eq. (2.23) is due to the matching 
with the preceding matter era. 
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